Zeno physics in ultrastrong circuit QED 
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We study the Zeno and anti-Zeno effects in a superconducting qubit interacting strongly and 
ultrastrongly with a microwave resonator. Using a model of a frequently measured two-level system 
interacting with a quantized mode, we predict different behaviors and total control of the Zeno times 
depending on whether the rotating-wave approximation can be applied in the Jaynes-Cummings 
model, or not. We exemplify showing the dependence of our results with the properties of the initial 
field states. 



I. INTRODUCTION 

Frequent measurements of a quantum mechanical sys- 
tem may cause its evolution to change, slow down or even 
freeze. This phenomenon, known as the quantum Zeno 
effect [l], 0|, implies in particular that an atom that is 
being periodically or continuously monitored cannot de- 
cay, or will do it at a slower rate. This suppression of 
spontaneous emission was first observed in trapped ion 
experiments Q , followed by observations in the suppres- 
sion of Rabi oscillations in cavity QED [3] , the decay of 
ultracold atoms [5|], and has been studied associated to 
photodetection in circuit QED 

Far from being just an interesting consequence of quan- 
tum mechanics, Zeno physics is also a useful tool in quan- 
tum control: it can be used to protect quantum infor- 
mation in certain subspaces 0, to suppress decoherence 
during quantum gates H, or even as a means of doing ef- 
ficient quantum search [9| . Interestingly, it has also been 
shown that frequent measurements can produce the op- 
posite effect [10, 11], enhancing the decay of a quantum 
system, in what is known as anti-Zeno effect. This other 
possibility has potential applications for cooling [l^] and 
has been experimentally demonstrated [13 1. 

Superconducting circuits are good candidates for 
studying Zeno physics [l^ . The equivalent of an atom 
interacting with the electromagnetic field in free space 
or in a cavity is replaced by a superconducting qubit 
interacting with an open or closed microwave transmis- 
sion line |l5|. By monitoring the state of the qubit, or 



of the cavity, one expects to have the possibility of freez- 
ing or accelerating the evolution. To our advantage, the 
time scales involved in the circuit QED dynamics and 
its measurements are now much larger as compared to 
the optical domain, and can be resolved with ordinary 
electronics. Furthermore, all energy scales of related ex- 
periments, including couplings, qubit and resonator fre- 
quencies, can be engineered and tuned at will. This 
opens the possibility of reaching the ultrastrong coupling 
regime |16l - [l9| . where the rotating-wave approximation 
(RWA) breaks down and the system cannot be described 
by a Jaynes-Cummings dynamics (20j . 

The goal of this work is indeed to study the Zeno 



physics in the qubit-cavity coupling of circuit QED with 
and without the RWA, that is, in the strong and ultra- 
strong coupling regimes. We will do it by using the 
model of a two-level system of frequency wq interact- 
ing with a single-mode harmonic oscillator of frequency 
Lo, and assuming periodic measurements of the qubit 
with intervals 5t. As we will show, there are two pos- 
sible regimes. First, the case where the measurement 
is faster than the vacuum Rabi coupling g of the com- 
bined qubit-resonator system, but slower than the qubit 
and resonator frequencies: g <C {5t)~^ ^ {u},ioq\. Here, 
the RWA is valid and we recover the usual slow down 
of the qubit decay. Second, the case where the qubit 
is measured more rapidly than any existing energy scale, 
{5t)^^ ^ {g, Lu, luq). Here, the RWA breaks down and we 
can observe both an enhancement of the Zeno or even an 
anti-Zeno effect, in which the decay of the excited state 
is accelerated with respect to the RWA case. 

The structure of the paper is as follows. We introduce 
the model of a qubit interacting with a single-mode res- 
onator and derive the evolution of the system under re- 
peated measurements. It will be shown that fast enough 
measurements induce an evolution of the qubit popula- 
tion that is well approximated by an exponential decay, a 
trace of the Zeno physics. The rate of the effective decay 
is shown to be different for a model that follows the RWA 
and for one that does not. We will discuss the conditions 
for these differences and explore their consequences using 
coherent and squeezed states, both analytically and nu- 
merically. Finally, we will suggest realistic implementa- 
tions of the proposed ideas in current setups of quantum 
circuit technologies. 



II. SURVIVAL PROBABILITY AND ZENO 
EFFECT 



Let us consider the dynamics of a two level system, 
g) and |e) , separated by a gap cjo, and coupled to a har- 
monic oscillator of frequency uj. The whole system models 
a superconducting circuit coupled to a high-Q microwave 
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resonator, and corresponds to the Hamiltonian {h = 1) 

H — gr (c+a + cr_at) + g^, (cr+a^ + (T_a) (1) 
t '^O 

+ uja'a+ -^<^z- 

The couphngs associated with the co-rotating Jaynes- 
Cummings (denote by "r" , from red) and counter- 
rotating Anti- Jaynes-Cummings (denote by "b" , from 
blue) are equal, g — gr — gh, but we express them sepa- 
rately to trace their physical consequences. It is known 
that when uj and Wq are much larger than the coupling 
strength, g, the RWA applies and formally g^ = 0, yield- 
ing the Jaynes-Cummings model ^20j] . 

We will focus on the survival probability, Pc{t): the 
probability of preparing the two-level system in a given 
state, say |e), and finding it unaltered at a given time t. 
Consider a system described by a qubit in an excited state 
and the cavity in an arbitrary state p/, that is p{0) = 
e) (e| eg) p/. The survival probability of this excited state 
at the first measurement is a function of time 



Pe(i)=Tr[p(i)le)(e|: 



(2) 



= Tr[e-*^V(0)e"^*|e)(e|] 

.»i=0 

where [H, p(0)]„ := [H, [H,... [H, p{0)]]] comes from the 
Baker-Campbell-Hausdorff formula. 

The short-time dynamics can be determined by the 
lowest order non-vanishing contribution in the previous 
sum. It can be shown that the linear term, 0{t), vanishes 
for arbitrary initial atomic states. Therefore the survival 
probability for short times decays at most quadratically 



(3) 



with a scale tz known as the "Zeno time" j21i] . If instead 
of performing just one measurement, one asks for the re- 
peated survival probability, Pe^^^NSt), where N measure- 
ments at regular intervals St were made, the result is no 
longer a quadratic but rather the exponential decay 



N 



-7cfft 



(4) 



The effective decay rate, 7ofr, and resulting lifetime. 
Toff, depend on the Zeno time 



7cff = 1/reff 



St/rl 



(5) 



and for frequent enough measurements, we will observe 
an effective freezing of the dynamics of the measured sys- 
tem, the proper Quantum Zeno effect 

The previous general considerations can be particular- 
ized for the Hamiltonian (IT|) to obtain the Zeno time 



g^gh {a- 
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FIG. 1. (Color online) Zeno effect in the strong coupling 
regime, g/u) — 0.05, for an initial state with an excited qubit 
and the resonator in a coherent state, |e) ® |a) , with lap — 9 
photons on average. In solid line we plot the collapses and 
revivals under free evolution in the RWA. The remaining lines 
correspond to the survival probability of the excited state un- 
der periodic measurements of the qubit, repeated at intervals 
St — O.OlTrg"^. If the RWA is applied, the usual Zeno time 
is obtained (dashed, red). If the RWA is not applied, we can 
either slow down the decay, setting the phase of the coherent 
state 5 = 7r/2 (green, dash-dot) or increase it 6 — (dotted, 
blue). For sufficiently fast measurements, these numerical re- 
sults and those obtained from the analytical expressions ([9]) 
and (jlOp are practically the same. 



where the expectation values {X) ~ Trt [p/AT] are 
taken over the initial state pf. Let us note that the Zeno 
time and the resulting decay rate are not only very sen- 
sitive to the initial field states, pf, but also to the precise 
form of the Hamiltonian. If the RWA is valid, g^ = 0, we 
obtain that the time depends on the total number of pho- 
tons, Tz^'^ oc (fltfl) + 1, while in the ultrastrong coupling 
regime it also involves other field quadrature momenta. 

The previous derivations are valid only under certain 
restrictions. In particular, the interval between measure- 
ments has to be short, so that we can replace the polyno- 
mial expression in @ with an exponential. This can be 
done if the time 5t is much shorter than the relevant pe- 
riods. In the case of a strong coupling regime with RWA, 
we will only need that this is shorter than the Rabi fre- 
quency, that is about St <^ l/ji. However, for ultrastrong 
couplings in which the RWA does not apply, we must 
impose faster measurements St <^ {l/"^, V.?}- 

It is also important to note that the Zeno time ([5]), 
and thus the effective decay rate ([5]), depends only on 
the initial field state pf, which is assumed constant. We 
may confirm this by performing a similar calculation as 
before. After N observations of the internal state of the 
two-level system every St, the field state will be 
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where Tr^ denotes the trace over the two-level system. By 
measuring the internal state of the two-level system and 
restricting ourselves to the survival probability, the field 
state is also frozen, i.e., the initial field state after each 
measurement is the initial field state plus some higher 
order corrections, as confirmed by the numerics. 

III. DIFFERENT INITIAL STATES 

Let us now study the changes in the Zeno effect de- 
pending on the initial resonator state, as well as the ex- 
tent of the differences between RWA and non-RWA mod- 
els. We begin with the Fock state, which has a well de- 
fined number of photons, \n) and gives an effective decay 
rate 

{l.s)r, = 5t[g^,(n + l)+gln] (7) 

that gives a difference of about a factor of two between 
the RWA and the non-RWA results. 

Much more interesting is the case when the field is 
assumed to be in a coherent state |a) [H, with mean 
photon number n and a phase 6*, that is a = n^/^e*^. The 
effective decay rate becomes 

(7eff)„='5i[.9?(l + l«n+5b|aP 

+ 2gr5b|a|'cos20] . (8) 

Out of the two experimentally accessible regimes we then 
have that 

(7eff)„ ='5t5'(l+4|a|'cos2(?), (9) 

{lfr'')^=5t9\l + \a\^). (10) 

Note that only when the RWA is not applied, the Zeno ef- 
fect becomes sensitive to the phase of the coherent state. 

Finally, we will consider a field initially in a vacuum 
squeezed state, l^/ij) , that is 

1^^) = e^^^*'''-^'^^') |0) , (11) 

where ^ — re'^ is the squeezing parameter. The decay 
rates, beyond and within the RWA, now become 

{'^cs)^ = 5t (cosh 2r - cos 6 sinh 2r) (12) 
(7cT^), = '5<.g2(coshr)^ (13) 

We have found that for coherent and squeezed states, 
the decay rate depends on the field phase. This depen- 
dence is enhanced in the non-RWA regime, which allows 
us for a total control of the decay rate. In the case of 
coherent states, we can make the decay rate independent 
of the number of photons. In the case of squeezed states 
we can change the way the decay rate depends on the 
squeezing, r: for the RWA or for = 7r/2 the Zeno decay 
increases as 1 + while for other values of 9 we may 
achieve 7cff = or make it decrease with the squeezing 
parameter as 1 — 2r if 61 = 0. We may regard the later as 
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FIG. 2. (Color online) Zeno effect in the ultrastrong cou- 
pling regime, when g, u), u)(i = 1 GHz. (top) Free evolution 
of the excited state population in the qubit and (bottom) 
survival probability of the excited state for periodic measure- 
ments with &t = Q.lu]~^ — 0.1ns (each symbol represents a 
measurement), (inset) Comparison between both plots. The 
solid line corresponds to a RWA model (black, solid, circles), 
while the other lines start from a coherent state with 9 pho- 
tons and either 6 = Q (blue, dashed, squares) or 6 = 7r/2 (red, 
dash-dot, triangles). 



an example of the quantum Zeno effect and the former 
as an example of the so-called anti-Zeno effect, where the 
dynamics of the system is accelerated rather than slowed 
down by frequent measurements. Note that this is an 
acceleration with respect to the RWA result, not to the 
natural decay rate since the system is not unstable. The 
transition from Zeno to anti-Zeno regimes in truly unsta- 
ble systems is discussed in Ref. [2^ . A possible definition 
of the Zeno and anti-Zeno time scales in oscillatory sys- 
tems (not unstable) is given in Ref. [22j |. 

The previous results can be numerically verified both 
in the strong and ultrastrong coupling regimes. In Figs.[T] 
we plot the evolution of a qubit-resonator system in 
which the coupling strength is small enough to satisfy the 
RWA in the free-evolution case, g = 0.05a; with uj = wq. 
Beginning with a coherent state in the field and an ex- 
cited state in the qubit, the results are collapses and re- 
vivals of the excited state population, a phenomenon that 
can be described with the Hamiltonian ([T|) with g\, — 0. 
However, if we perform measurements with a periodicity 
5t — 0.27rcj^^ ~ O.OItt^^^, we observe that the survival 
probability of the excited state follows an exponential 
law ([4]) with a decay rate that is not the one given by the 
RWA (red, dashed line), but rather the one in Eq. ©. In 
particular, the decay rate can be enhanced by setting the 
resonator in a coherent state with = 0, or decreased, 
by choosing 6 = 7r/2. A similar result is observed for 
ultrastrong couplings, g = uj = luq, as in Fig. [5J 
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FIG. 3. A phase qubit can be used to explore the Zeno effect: 
the escape rate of the unstable level |e) is equivalent to a 
continuous measurement of that state. 

IV. EXPERIMENTAL IMPLEMENTATION 

The main question that arises when looking for an ex- 
perimental implementation of the previous ideas is how 
to achieve the regime of fast repeated measurements in 
superconducting circuits. As discussed before, we need 
to probe both a regime in which the time interval St 
is much smaller than the vacuum Rabi period but the 
RWA remains valid, and a regime in which the measure- 
ments take place much faster than any other time scale 
of our system. Given that current measurement devices, 
as SQUIDs, need nanoseconds to operate, we would have 
to relax all energy scales, oj, ujq, and and risk being 
affected by thermal fluctuations. However, this is not un- 
conceivable because we could design an experiment near 
the largest frequency w = 2tt x 200 MHz, which gives a 
vacuum Rabi period of about 5 ns. The cavity will then 
be populated by an average of 2 thermal photons, but 
the features of the Zeno effect should remain visible. 

Another possibility is to work in a framework of con- 
tinuous measurements, like in fluorescence experiments 
with trapped ions fsj. In that case, we suggest the use 
of a resonator interacting with a phase qubit [28| whose 
excited state, |e) may decay into the continuum at a rate 
r, [See Fig. [3]. By studying the repeated survival proba- 
bility of the qubit ground state, we should observe also a 
Zeno effect, where the incoherent decay plays the role of 
a continuous measurement. According to Refs. j26l. [27|. 



the decay channel of the qubit excited level gives an ef- 
fective measurement interval St = A/T. By continuously 
monitoring whether the phase qubit switches to a voltage 
state, we will be able to determine the survival probabil- 
ity of the ground state. The phenomenology should be 
equivalent to the previous cases, with the advantage that 
now Lu can be larger, of the order of GHz, due to the 
flexibility in tuning the decay rate T. 

Arguably, the most promising setup would consist on a 
superconducting qubit interac ting via a tunable coupling 
with the microwave resonator |29| . By controlling the in- 
teraction we should be able not only to choose between 
strong and ultrastrong couplings, but also to completely 
deactivate the interaction for brief periods of time in or- 
der to measure the system. This switching on and off can 
be done with sub-nanosecond resolution, allowing us to 
reach the regime St ^ {1/cj,1/cjo,1/5}, and the whole 
system should behave just as if we were able to instanta- 
neously monitor the superconducting qubit. 

Summing up, in this work we have shown the difference 
in the Zeno dynamics between a quantum optical system 
that obeys the rotating- wave approximation and one that 
does not. In the RWA case, the Zeno effect only depends 
on the total number of photons of the electromagnetic 
field that interacts with the monitored qubit. In the non- 
RWA case, the Zeno decay may be completely suppressed 
or enhanced (anti-Zeno), by suitably preparing the phase 
of a coherent or squeezed field. Both regimes should be 
observable in quantum circuits consisting on a qubit in- 
teracting with a high-quality microwave resonator. Note 
that in contrast to previous works (sO] , the choice of the 
physical system makes the discrimination of RWA/non- 
RWA more accessible and also avoids the subtleties and 
problems associated with a continuum of modes. 

The authors acknowledge funding from UPV-EHU 
Grant GIU07/40, Spanish MICINN project FIS2009- 
12773-C02-01, EuroSQIP and SOLID European projects, 
ESF network program "AQDJJ", Basque Government 
Grants BF108.211 and IT472-10, Spanish MEG Project 
FIS2006-04885 and CSIC Project 2008501044. 



[1] B. Misra and E. Sudarshan, J. Math. Phys. 18, 756 
(1977). 

[2] P. Facchi, S. Pascazio, J. Phys. A 41, 493001 (2008). 
[3] W. M. Itano, D. J. Heinzen, J. J. BoUinger, and D. J. 

Wineland, Phys. Rev. A 41, 2295 (1990). 
[4] J. Bernu, S. Deleglise, C. Sayrin, S. Kuhr, I. Dotsenko, 

M. Brune, J. M. Raimond, and S. Haroche, Phys. Rev. 

Lett. 101, 18 (2008). 
[5] E. W. Streed, J. Mun, M. Boyd, G. K. Campbell, P. 

Medley, W. Ketterle, and D. E. Pritchard, Phys. Rev. 

Lett. 97, 260402 (2006). 
[6] F. Helmer, M. Mariantoni, E. Solano, and F. Marquardt, 

Phys. Rev. A 79, 052115 (2009). 
[7] A. Barenco, A. E. Berthiaume, D. Deutsch, A. Ekert, R. 



Jozsa and C. Macchiavello, SIAM J. Comput 26, 1541 
(1997). 

[8] A. Beige, D. Braun, B. Tregenna, and P.L Knight, Phys. 

Rev. Lett. 85, 1762 (2000). 
[9] A. M. Childs, E. Deotto, E. Farhi, J. Goldstone, S. Gut- 

mann, and A. J. Landahl, Phys. Rev. A 66, 032314 

(2002). 

[10] A. G. Kofman, and G. Kurizki, Nature 405, 546 (2000). 
[11] K. R. Chapin, and M. O. Scully, Nature 452, 705 (2008). 
[12] N. Eretz, G. Gordon, M. Nest, and G. Kurizki, Nature 

452, 724 (2008). 
[13] M. C. Fischer, B. Gutierrez-Medina, and M. G. Raizen, 

Phys. Rev. Lett. 87 , 040402 (2001). 
[14] J. Gambetta, A. Blais, M. Boissonneault, A. A. Houck, 



5 



D. I. Schuster, and S. M. Girvin, Phys. Rev. A 77, 012112 
(2008). 

[15] A. Wallraff et al, Nature 431, 162 (2004). 

[16] J. Bourassa, J. M. Gambetta, A. A. Abdumalikov Jr., O. 

Astafiev, Y. Nakamura, and A. Blais, Phys. Rev. A 80, 

032109 (2009). 

[17] A. A. Abdumahkov, Jr., O. Astafiev, Y. Nakamura, Y. A. 

Pashkin, and J. S. Tsai, Phys. Rev. B 78, 180502 (2008). 
[18] T. Niemczyk, F. Deppe, H. Huebl, E. P. Menzel, F. 

Hocke, M. J. Schwarz, J. J. Garci'a-Ripoll, D. Zueco, T. 

Hmmer, E. Solano, A. Marx, R. Gross. .,arXiv:100 3.2376 
[19] P. Forn-Di'az, J. Lisenfeld, D. Marcos, J. J. Garcia- 

Ripoll, E. Solano, C. J. P. M. Harmans, and J. E. Mooij, 

arXiv:1005.1559 
[20] E. T. Jaynes and F. W. Cummings, Proc. IEEE 51, 89 

(1963). 

[21] L. Schulman, A. Ranfagni, and D. Mugnai, Phys. Scr. 



49, 536 (1994). 

[22] P. Facchi and S. Pascazio, Quantum Zeno and inverse 
quantum Zeno effects, Progress in Optics, Vol. 42, Chap- 
ter 3, edited by E. Wolf (Elsevier, Amsterdam, 2001). 

[23] D. F. Mundarain and J. Stephany, Phys. Rev. A 73, 
042113 (2006). 

[24] D. Mundarain, M. Orszag, and J. Stephany, Phys. Rev. 

A 74, 052107 (2006). 
[25] K. Modi and A. Shaji, Phys. Lett. A 368, 215 (2007) 
[26] L. S. Schulman, Phys. Rev. A 57, 1509 (1998). 
[27] J. G. Muga, J. Echanobe, A. del Campo, and I. Lizuain, 

J. Phys. B 41, 175501 (2008). 
[28] J. M. Martinis, M. H. Devoret, and J. Clarke, Phys. Rev. 

Lett. 55, 1543 (1985). 
[29] B. Peropadre, E. Solano, and J. J. Garcia-Ripoll, 

arXiv:0912.3456 
[30] H. Zheng, S. Y. Zhu, and M. S. Zubairy, Phys. Rev. Lett. 

101, 200404 (2008). 



